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Abstract 


Using a temperature-dependent potential obtained from lattice gauge cal¬ 
culations of Karsch et al., we study the stability of heavy quarkonia in the 
quark-gluon plasma. We find that only the T(IS') and the ^(IS) are bound 
in the quark-gluon plasma, and have a small binding energy. The quark- 
gluon plasma may be revealed by an T(IS') dilepton peak with an invariant 
mass close to twice the current b quark mass, which is lower than the T(IS') 
mass in free space. The quarkonia T(IS') and r]b{lS) can dissociate by colli¬ 
sion with quarks and gluons in the quark-gluon plasma. The T(IS') and the 
r]b{^S) can also dissociate spontaneously at temperatures above the dissocia¬ 
tion temperature l.llTc, where Tc is the quark-gluon plasma phase transition 
temperature. At temperatures slightly above the dissociation temperature 
these states appear as resonances, which provides another signature for the 
quark-gluon plasma. 
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I. INTRODUCTION 


The suppression of heavy quarkonium production in a quark-gluon plasma has been a 
subject of intense interest since the pioneering work of Matsui and Satz |l[]. Initial insight 
into the dissociation temperatures of heavy quarkonium was further provided by Karsch, 
Mehr, and Satz [0. Recently, using the temperature-dependent Q-Q interaction inferred 
from lattice gauge calculations of Karsch et al. [^, Digal, Petreczky, and Satz reported 
theoretical results for the dissociation temperatures of heavy quarkonia in hadronic matter 
(T < Tc) and in the quark-gluon plasma (T > Tc), where Tc is the quark-gluon plasma 
phase transition temperature. Subsequent analysis of the dissociation temperatures of heavy 
quarkonia at T < Tc using different interactions and selection rules was given in detail in 
1^ and summarized in [^. The dissociation of heavy quarkonia by thermalization and by 
collisions with particles in the medium has also been investigated 

For T > Tc, the equilibrium state of matter is the quark-gluon plasma. If the interaction 
between a QQ and the QGP is weak, the QQ system can be color-singlet or color-octet. 
Using perturbation theory as a guide, Digal et al. extracted the color-singlet potential 
Vi{r,T) and color-octet potential Vs{r,T) of a heavy quark pair from the Polyakov loop 
(T(O)T'^(r)). They used the relation 
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where (T^) denotes (T(O)T^(r)) at r —>■ oo. This relation is valid when the QQ interacts 
weakly with a perturbative color medium so that color-singlet and color-octet QQ states are 
approximate eigenstates of the system. They also made the assumptions 
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as suggested by perturbation theory (for which C(T) = 1). They assumed that /i(T) was 
1.15T from the systematics at high temperatures and chose functions q:{T) and c(T) to £t 
the Polyakov loop results of Karsch et al. [^]. They examined color-singlet heavy quarkonia 
in the quark-gluon plasma, even for temperatures slightly greater than Tc 0. 

We shall restrict our interest to the region of temperatures slightly greater than Tc, which 
is the most important region for the investigation of the stability of heavy quarkonia in the 
quark-gluon plasma. In this temperature region, the screening effects are non-perturbative. 
Evidence for non-perturbative screening comes from numerical lattice gauge calculations 
of Karsch et al. 0 who show that the Debye screening mass is significantly larger than 


the leading-order perturbative value of = gT^l + nf /6 and that the screening mass 
remains about three times larger than the perturbative Debye mass even up to temper¬ 
atures as high as lOOTc. Similarly, the rT dependence as shown in Eqs. 0 - 0 ) with the 
perturbative value of C(T) = 1 is approximately valid only for T > 3Tc pl| . 
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For the region of temperatures close to Tc in which we shall be interested, C{T) obtained 
from the procedures of Ref. is about 0.1 which deviates signihcantly from the perturbative 
value of C{T) = 1, and the system cannot be described as a weakly coupled and pertur- 
batively screened system [^. Non-perturbative screening effects are important and the 
extrapolation procedure in the non-perturbative region employed in Ref. is ambiguous. 
Specihcally, the Polyakov loop provides a single number for a given r and T, but the above 
separation procedure yields quantitative determination for the physical quantities of Ri, Vs, 
and the (dynamical) relative weight of the color-singlet and color-octet components. Many 
assumptions are needed to extract these three physical quantities out of a single quantity 
of the Polyakov loop. Without theoretical non-perturbative results to guide us, these ad¬ 
ditional assumptions are not unique. Different assumptions will lead to different extracted 
quantities in the non-perturbative region. 

Related to the question of the separation of the potential is the question whether purely 
color-singlet and color-octet QQ states are approximate eigenstates in the non-perturbative 
region of temperatures in QGP, when T is just slightly greater than T^. Purely color-singlet 
and color-octet QQ states are approximate eigenstates of the system in the high-temperature, 
perturbative region for which Eq. ([1.1|) is valid. This however cannot be the case in the non- 
perturbative region of T> Tc- A constituent of a QQ in a deconhned medium can interact 
with a constituent of the medium (g, q, and g) by exchanging a single gluon. Such a gluon 
exchange alters the color of the QQ from a color-singlet state to a color-octet state and 
vice versa. The color of the medium will undergo corresponding changes so that the color 
of the whole system is neutral. Such a color-changing interaction mixes the color states 
of the heavy quarkonium. We shall show below in Section II that the degree of mixing 
increases with the magnitude of the ratio of the color-changing interaction matrix element 
A to the energy difference A between the unperturbed color-octet and color-singlet states. 
In the non-perturbative region of T in the QGP, |A/A| is not small and the eigenstates of 
the system become states of mixed colors. A purely color-singlet or purely color-octet state 
cannot be a stationary state of the non-perturbative system. 

Even though a QQ state in the non-perturbative region of temperatures cannot be sep¬ 
arated into eigenstates of pure color, it is still a meaningful question to hnd out how the Q 
interacts with the Q in such a system, whatever the color admixture of the resultant QQ 
may be. In lattice gauge calculations the expectation value of the Polyakov loop operator 
is evaluated allowing all possible configurations of the gauge helds and fermion fields at all 
lattice points and links, and this corresponds to all color states of the QQ and the medium. 
The gauge held variables along the path of the Polyakov loop provide information on the 
evolution of the heavy-quark system. The greatest contributions to the Polyakov loop at 
a given Q-Q separation will come from those gauge helds and quark helds which adjust 
themselves in magnitude and in color admixture to give states of lowest energies of a Q-Q in 
the medium. As the eigenstates of the system are states with mixed colors (see Section II), 
the Polyakov loop provides information on the free energy of mixed-color states. Garrying 
this calculations for diherent separations between the Q and the Q, one obtains the free 
energy as a function of the separation, which provides the spatial variation of the ehective 
interaction between the Q and the Q. It can be used to hnd out if the interaction is strong 
enough to bind the Q with the Q in the non-perturbative deconhned medium. Therefore, 
in this non-perturbative medium with strong interactions, the ehective interaction potential 
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VQQ(r, T) between the Q and the Q, is related to the expectation value of the Polyakov loop 
from lattice gauge theory by 


{L^> 


= exp 


-VQQir,T) 

T 


(1.4) 


The quantity VQQ{r,T) was extracted previously from the lattice gauge calculations by 
Karsch et al. 0. The binding energies of quarkonia as a function of the temperature will 
provide useful plasma diagnostic information on the dissociation temperatures and quarkonia 
masses in the nonperturbative quark-gluon plasma phase. 

As a heavy quarkonium becomes unbound, it appears as a resonance in the continuum. 
Such a resonance occurs only in the quark-gluon plasma at a certain range of temperatures 
and may be used as a signature of the quark-gluon plasma Following the method used 
previously by Calogero and by Wong and Chatterjee |^, we shall study the resonance 
structure of heavy quarkonia in the quark-gluon plasma. 

Quarkonia can also dissociate by collision with particles in the medium In par¬ 

ticular, we will study T(IS') + g —>• b + b. Such a dissociation process is similar to the 
photodissociation of a deuteron after absorbing a photon. This dissociation process can be 
studied by using a multipole expansion of the gluon held, analogous to the electromagnetic 
case discussed in detail by Blatt and Weisskopf |]2^. As the gluon dissociation cross section 
has also been obtained previously by Peskin and Bhanot p5| , ^ and subsequently used by 
Kharzeev and Satz |]T^, we show in the Appendix the equality of the cross sections obtained 
by the two formalisms for the simple case considered by Peskin and Bhanot. 

In Section II, we discuss the color structure of a QQ and its dependence on the interaction 
between the constituents of the QQ and the constituents of the medium. In Section III, 
we describe the Schrodinger equation used to calculate the energies and wave functions of 
quarkonium states. We introduce an analytical form of the effective screening potential 
Vqq to represent the results of lattice gauge calculations. The quarkonium energies and the 
dissociation temperatures are calculated in Section IV. In Section V, we consider continuum 
states and follow the resonances as functions of the quark-gluon plasma temperature. In 
Section VI we evaluate the cross section of the dissociation process T(IS') + g ^ b + b. A 
discussion of our results is given in Section VII. 


II. COLOR STRUCTURE OF A QQ IN THE QGP 

The color structure of a QQ placed in a deconhned medium depends on the interaction 
between the constituents of the QQ and the constituents of the deconhned medium. We can 
consider hrst two extreme cases of very weak and very strong interactions as illustrated in 
Fig. 1. 

If the interaction between the Q and the Q of the quarkonium with q, q, and g is weak as 
in a perturbative medium, [Fig. 1(a)], the Q and the Q can combined together to form pure 
color-singlet and color-octet eigenstates. The concept of a purely color-singlet or color-octet 
state is useful as an approximate description, when the interaction between the Q and the 
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Q with the quark matter particles can be ignored. 


(a) Perturbative Medium 
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(b) Highly Non-perturbative Medium 



Fig. 1. Q and Q in a medium with perturbatively weak interactions 
(a), and in a highly non-perturbative medium with strong interac¬ 
tions (b). Open circles represent quarks, antiquarks, and gluons in 
the plasma. 

In the other extreme of very strong interactions as illustrated in Fig. 1(b), the Q and the 
Q are intimately linked to the medinm particles by gluons. The isolation of the QQ system 
is impossible and one cannot speak nnambiguously of a purely color-singlet or color-octet 
state. 

We discuss how we can pass from one limit to the other limit with a schematic model. 
We consider a QQ system in a medinm and use {QQY^^ basis states where c = 0 is for 
color-singlet and c = 8 for color octet. We consider for simplicity the medinm represented 
by where c is again the color index. As the color of the whole system consisting of both 
{QQY and must be color neutral, a state of the system is described by 

1 ^) = + as[{QQ)^^^M^^Y^ ( 2 . 1 ) 

One can represent this admixture by a column vector with elements oq and og. In the 
absence of the interaction between the Q and the Q with the medinm particles, the color- 
singlet state and the color-octet state can be separated and we denote the energy difference 
between the unperturbed color-octet state and the color-singlet state by A = — E^^'>. 
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By the interaction VAai^ + o, A' + a') in which the constituent A in exchange 

a gluon g with the constituent a in M^) to become A' and a', the color of both and 

will change. We define the color-changing transition matrix element A as 

A = ([(gg)(«)M®]W \vAa\ (gg)(°)M(°)) (2.2) 


The Hamiltonian for such a system is 

_ / -A/2 A 
A A/2 

The lowest energy eigenstate of the system is then 
A + (1 + 4AVA)2)'/2'j 
= [ -2A/A 

The excited eigenstate is 

2A/A 


jv^|l 


+ 4AVA2 + (1 + 4AVA2)i/2]i/2 
1 


= 


+ (1 + AxyAyy/y ^/2[l + 4A7A2 + (1 + 4A2/a2)i/2]i/2 


(2.3) 


(2.4) 


(2.5) 


Thus, the color-structure of the system depends on the ratio |A/A|. In the confined 
phase, A is very large, and so |A/A| is very small; the lowest eigenstate is purely color 
singlet. In the perturbative deconhned phase at very high temperatures, A is small and 
I A/AI is also small; the eigenstates are nearly pure in color. 

On the other hand, if |A/A| is large, the eigenstates of the system is an admixture of 
color-singlet and octet state. The nonperturbative region of temperature is a case of strong 
interactions between the Q (or Q) and the medium particles. We can estimate the ratio 
I A/AI which determines the color admixture. The energy difference A between color-octet 
and color-singlet states is approximately the dissociation energy of the color-singlet state. 
Earlier estimates of the dissociation energy of J/ip in the hot hadron phase show that it 
decreases with increasing temperatures, resulting in zero dissociation energy or nearly 
zero dissociation energy (see Fig. 3 of @]) for J/tp at T< 7 . For T> 7, screening of 
deconhned quark matter will tend to weaken further the c-c attraction. (Calculations in 
Section IV is consistent with a vanishing value of the dissociation energy of J/pJ.) Thus, A 
is zero or nearly so for cc. For bb, the dissociation energy is about 0.15 GeV for T(IS') 0 
at T7 Tc- For T7 7, screening in the quark-gluon plasma makes the dissociation energy 
even smaller. (In Section VI, the dissociation energy is found to be of order 0.01 GeV.) 
On the other hand, for a QGP at T = 200 MeV, the average separation between the 
quanta in the QGP is approximately d= 0.6 fm . The transition matrix element is of order 
A ~ exp{—mthd}/d where the thermal mass rrith ~ gT and g"^ = as/4:7i. We hnd that 
for as = 0.2, A is of order 0.06 GeV. The ratio |A/A| is also quite large for bb, if we take 
A ~ 0.01 GeV. Therefore, the QQ system is strongly mixed in color in the nonperturbative 
deconhned region of temperatures slightly greater than Tc. 

Even though a gg state in the non-perturbative region of temperatures cannot be sep¬ 
arated into eigenstates of pure color, it remains a meaningful question to hnd out how the 
g interacts with the Q in such a system. As explained in the Introduction, the expectation 
value of the Polyakov loop operator gives the free energy of the system of the QQ coupled to 
the medium in a mixed-color state. It provides the ehective interaction potential VQQ{r,T) 
between the Q and the Q in in the region of T> T^. We shall use this interaction potential 
to examine the stability of the gg system in the next few sections. 
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III. SCHRODINGER EQUATION FOR HEAVY QUARKONIUM STATES 


The energy e(T) of a heavy quarkonium state {QQ)jls, measured relative to the mass of 
the Q and the Q, can be obtained by solving for the associated eigenvalue of the Schrodinger 
equation 


+ VQQ{r,T)^^l:jLs{r,T) = emjLs{r,T), (3.1) 

where fiqq = mq/2 is the reduced mass. The quarkonium is bound if e(T) is negative, and 
is unbound and dissociates spontaneously into a Q and a Q if ^(T) is positive. 

The color charges of the Q and the Q are screened in the quark-gluon plasma. We can 
conveniently represent Vqq{r,T) by a screened Yukawa potential with an effective strength 
aeff(T) and a screening mass parameter /i(T). We shall be interested in temperatures slightly 
greater than Tc. The Q-Q interaction obtained from the lattice calculations of Karsch et al. 
0 for 1.2Tc > T > Tc in the quark-gluon plasma phase can be approximated by 

4 p-KT)r 

= --aefr(T)—-— (3.2) 

where 

«efr(7") = 0.20, /i(T) =/xq + hT(7"/Tc — 1), (3.3) 

Po = 0.323 GeV, and /xt = 3.04 GeV. 

Fits to Vqq{r,T) using the above form are shown in Fig. 2. We include the color-singlet 
factor (—4/3) in the above definition (|3.2|) for Q;efr(^) m order to have a comparison of the 
strength of the Q-Q interaction in the quark-gluon plasma phase with the color-singlet Q-Q 
interaction seen at T < Tc p|. The strength parameter aeff(T) has very weak temperature 
dependence, and can be adequately represented by the constant aefr(T) = 0.20. This value 
of aeff(T) = 0.20 in the quark-gluon plasma phase is smaller than the value of as ~ 0.32 
for charmonium and ag ~ 0.24 for bottomium in the screened color-Goulomb potential for 
temperatures below Tc [§. The screening mass parameter /xq for T > Tc is greater than 
the screening mass parameter of /xq = 0.28 GeV used for T < Tc in p. One concludes 
that because of the screening due to deconfined quarks and gluons, the effective potential 
between the Q and the Q in the quark-gluon plasma remains attractive, but its strength is 
weaker and its range is shorter than the color-Goulomb interaction seen in hadronic matter. 
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The range of the screened potential decreases with increasing temperatnres. 



Fig. 2. The effective potential VQQ{r,T) between the Q and the 
Q in a qnark-glnon plasma at various temperatures. The symbols 
represent the results of the lattice gauge calculations of |P] reported 
in [^. The curves in the main figure are from the paramterization 
of Eqs. (|3.2|) and ( p.3|) . The inserted hgure shows the effective 
potential at T = 0(for charmonium) and T = 1.03Tc. 


IV. BOUND STATES AND DISSOCIATION TEMPERATURE 


The eigenvalues of the Hamiltonian can be obtained by integrating the Schrodinger equa¬ 
tion numerically and requiring an exponentially decaying wave function at large distances. 
For our numerical calculations we use a current c quark mass rric = 1.25 GeV, and a current 
b quark mass mb = 4.2 GeV . 


We hnd that the potential of Eqs. ( |3.2| ) and ( p.3|) , as determined by the lattice gauge 
calculations of Karsch et al., is too shallow to hold a bound charmonium state. Because 
of the large b quark mass, we can neglect the spin-orbit and hyperfine interactions. In this 
approximation, the L = 0 T(IS') and rjbi^S) states are degenerate, as are the L = 1 Xbi^P) 
and hb{lP) states. For bb the L = 1 Xb{^P) and hb{lP) states are unbound, and the only 
bound states are the L = 0 T(IS') and ? 7 b(lS') states. The energies of the T(IS') and rib{lS) 
states as functions of the quark-gluon plasma temperature are shown in Fig. 3. They become 


















less bound with increasing temperature. The dissociation temperature of T(IS') and rib{lS) 
is Td = l.llTc, as determined by the condition e = 0. The root-mean-squared radius of 
T(IS') at T = Tc is 0.80 fm, indicating a large radius because of the weak binding. The 
rms radius increases with temperature. At T = l.lOTc which is close to the dissociation 
temperature, the rms radius has increased to 3.3 fm. 



Fig. 3. The energies of T(ls) and rib{lS) states as a function of 
temperature. 

Previously we found that the Xbi and Xb 2 states were stable against dissociation for 
T < Tc- In the quark-gluon plasma phase, however, screening by the deconhned q, q, and g 
alters the interaction between Q and Q. The interaction does not need to be a continuous 
function of the temperature at T^, because a new type of screening arises from deconhned 
q, q, and g in the quark-gluon plasma phase for T > Tc. Consequently, the quarkonium 
energies do not need to be continuous across the QGP phase transition boundary. 

For completeness we show the heavy quarkonium dissociation temperatures in Table I 
where we also include those shown in Table III of 0. The dissociation temperature of 
Td/Tc = 1.00-1- for Xh and hb in Table I denotes Td/Tc = 1.00 in the quark-gluon plasma 
phase. The dissociation temperatures obtained here are substantially lower than the values 
of 1.13Tc for the color-singlet Xb and 2.31Tc for the color-singlet T(IS') obtained by Digal 
et al. using the color-singlet potential Vi(r, T). These differences arise because the present 
investigation assumes a QQ pair that is immersed in and interacting strongly with the color 
medium of the quark-gluon plasma, and the work of Digal et al. p] considers the dissociation 
of a color-singlet QQ pair in the quark-gluon plasma. Hence there are significant differences 
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in the dissociation temperatures. 


Table I. The dissociation temperatures in units of for various heavy quarkonia 


Heavy 

Quarkonium 


Xc2 

Xcl 

J/^ 


Xb2 

Xw 

T'^ 

Xbi hb 

y^vb 

Td/n 

0.50 

0.91 

0.90 

0.99 

0.57 

0.82 

0.82 

0.96 

1.00+ 

1.11 

Td/% 

(Digal et al. [Q) 

0.1-0.2 

0.74 

0.74 

1.10 

0.75 

0.83 

0.83 

1.10 

1.13 

2.31 


V. QUARKONIUM STATES IN THE CONTINUUM 


To study the decay of a heavy quarkonium resonance, it is of interest to examine the 
related Q-Q scattering process. The scattering phase shifts with the screened potential 
VQQ{r,T) will provide information regarding the locations and the widths of heavy quarko¬ 
nium resonances in the continuum. 

To obtain the continuum wave function in the scattering problem, we use the phase-angle 
method discussed in detail by Calogero [^] and used previously by Wong and Chatterjee 
We write the wave function as 


‘ipjLsir) = RjLs{r)yjLs{r) = -^j^^yjLs{r) 


I dvr 
2 / + 1 ’ 


(5.1) 


and represent it in terms of an amplitude a^lr) and a phase shift S^lr) 


ujLsir) = 


adr) 

aiioo) 


Diihr) sm{6L{kr) + 6L{r)), 


(5.2) 


with the boundary condition that —>■ 0) = 0. We have normalized the radial wave 

function such that it coincides with the normalization of Eq. (4.32) of Blatt and Weisskopf 


2^. In the case of a free particle, 'Ujls'(^) ~ brj^ihr). The functions D{kr) and 6{kr) are 

(5.3) 


known functions determined from the free particle wave functions |2^ , 

IDq{x) = 1 , lDi{x) = (1 l/x^fR, 

and 

5 q{x) = x^ 5 i{x) = X — X. 


(6.4) 


The phase shift function Siir) depends on the interaction f/(r) = 2/rQQV'gg(r). The equation 
for ^^(r) is given by ^2 




f/(r) 


Dl{kr) |sin[(5L(/cr) -F (5L(r)]| 


(5.5) 


After the function Sii^r) is evaluated, the amplitude ai^r) can be obtained from ^^(r) by 


aL{r) = exp 


1 r 

2k Jo 


dsU{s) D\{ks) sin2[(5 l(/):s)-|-5l(s)] 


(5.6) 
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Fig. 4. (a) The phase shift (fg and (b) sin^ for fhe L = 0 

scattering of Q and Q as a fnnction of energy e for different T 
valnes near the dissociation temperatnre T^. 


Nnmerical integration of Eq. (|5.5|) gives the asymptotic phase shift 6l = 5l{oo) and the 
continnnm wave fnnction. We show the phase shift (5o in Fig. 4(a) and sin^ (5o in Fig. 4(6) 
for the L = 0 66 state as a fnnction of energy e at varions temperatnres. The phase shift 
6 o eqnals vr at e = 0 for T < 1.1 ITc, in accordance with the Levinson’s theorem [which 
states that 6 = rni at e = 0, where n is the nnmber of bonnd states, except when an S-wave 
resonance occnrs at zero energy, in which case 6 = {n + l/2)7r]. At temperatnres below 
the dissociation temperatnre = l.llTc the phase shift decreases as a fnnction of energy, 
and is eqnal to 7r/2 at some energy e. Snch an occnrrence of 60 = 7r/2 with d6{e)/de < 0 
does not represent a resonance; it is an “echo” of the bonnd state which lies jnst below the 
continnnm At the location sin^ 60 = 1 the scattering cross section is enhanced when a 


bonnd state lies jnst below the continnnm, as was hrst noted by Wigner and discnssed by 
Landan and Lifshitz . 

At the T(IS') dissociation temperatnre Td = l.llTc, 60 = 7 r /2 at e = 0, and the L = 0 
T(IS') and rjbi^S) resonances occnr at e = 0. The phase shift decreases fairly rapidly as a 
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function of energy. The energy difference between the location of the maximum of sin^ (5o 
and the half maximum is about 0.2 GeV, indicating a resonance with a half width of about 
this magnitude. 

For temperatures higher than the phase shift is zero at e = 0. It increases rapidly with 
e and reaches a maximum before it decreases slowly. One can associate the location e of the 
maximum phase shift in e with the position of the T(IS') and r]b{lS) “resonances”. Strictly 
speaking the enhancement at the phase shift maximum does not represent a resonance, as 
the magnitude of the maximum phase shift is less than 7r/2. However when the transition 
temperature is slightly greater than T^, the magnitude of the maximum phase shift is close 
to n/2, and it decreases continuously as the temperature increases. The location e of the 
maximum phase shift is a smooth continuation of the bound state energy. We can therefore 
speak of these enhancements as “resonances”, as they are the continuation of the bound 
states into the positive energy domain. We can choose to identify the position of the max¬ 
imum of the phase shift for T > with the location of the resonances as a function of 
temperature. The maximum value of sin^ (5o decreases as the temperature increases. The 
width of sin^ (5o increases, showing a broader enhancement of the cross section as the tem¬ 
perature increases. For example, at T = Td + 0.02Tc the maximum occurs at e = 0.008 GeV, 
and the separation between the maximum and the half maximum on the higher energy side 
is about 0.30 GeV. At T = + O.OdT^ the maximum occurs at e = 0.018 GeV, and the 

separation between the maximnm and the half maximnm on the greater energy side is about 
0.43 GeV. The function sin^ (5o is not symmetrical with respect to the maximum. It initially 
rises rapidly, but then decreases slowly as a function of e. The width increases rapidly with 
temperatnre. 

We show in Fig. 5 the phase shift 5i and sin^ for L = 1 Q-Q scattering. For T slightly 
greater than the phase transition temperatnre, the phase shift has a maximum at e ~ 0.1-0.2 
GeV. The maximum magnitude of the phase shift is 0.54 at e = 0.12 GeV for T/Tc = 1.01, 
and is 0.39 at e = 0.265 GeV for T/Tc = 1.05. The maximnm phase shift does not reach 
n/2. The qnantity sin^ (5i has a maximum value of 0.22 at T /Tc = 1.01, and a valne of 0.14 
at T/Tc = 1-05. These enhancements are qnite broad. For example, for T/Tc = 1-01, the 
qnantity sin^ 6i decreases slowly after it reaches the maximnm valne, and the separation 
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between the maximum and the half maximum on the higher energy side is about 1.7 GeV. 



£ (GeV) 


Fig. 5 . (a) The phase shift (5i and (b) sin^ (5i for the scattering of 
Q and Q in the L = 1 state as a function of the energy e for various 
T(T,. 


VI. CROSS SECTION FOR T(15) + b + b 

The quarkonia T(IS') and r]i,{lS) are stable in the quark-gluon plasma, but have small 
binding energies (from 0 to about 0.0175 GeV). They can dissociate into b and b by collision 
with quarks and gluons. 

We shall focus our attention on the dissociation reaction T(IS') +g b+b through gluon 
absorption; the dissociation of ^7(15') can be treated in a similar way. For the dissociation 
processes T(IS') + g ^ b + b, the initial quarkonium state is a IS* bound state, and the hnal 
state is a bb state in the continuum. This process is similar to the photodissociation of a 
deuteron after photon absorption. The dissociation process can be studied using a multipole 
expansion of the electromagnetic or chromodynamic helds, and the lowest multipole orders 
are Ml and El absorption processes. For low energy photons or gluons, these processes can 
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be treated in the long wavelength limit, discussed in detail by Blatt and Weisskopf . 
Following Eq. (XIL4.33) of Blatt and Weisskopf 
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the El dissociation cross section 
for a transition from the initial bound IS* state with energy e, = —B = to a hnal 


IP continuum state with energy e = k'^/mg, after absorbing a gluon of energy e- 
by 


'^dis 


= 4x|agQ(A;2 + 7^)A;-l/^ 


■Cj, IS given 


( 6 . 1 ) 


where the factor of 4 on the right hand side has been included for color electric dipole 
contributions from both the Q and the Q, (as distinct from the deuteron case, in which only 
the proton contributes to the electric dipole matrix element). 

The quantity I of Eq. O is the radial overlap integral for the El electric dipole matrix 
element. 


/ = / Uip{r)uis{r)rdr. ( 6 . 2 ) 

Jo 

where the wave function Mip(r) is given by Eqs. (|^) - (|5.6|) . The IS* bound state wave 
function in Eq. ( |6.2|) is normalized according to 

poo 

/ \uis{r)\^dr = 1. (6.3) 

Jo 

The quantity OgQ is an effective strong interaction coupling constant for a gluon interacting 
with a heavy quark or antiquark. Specihcally, similar to Eq. (^ of the Appendix, we have 

dvragQ = |(/|^y|i)|^ (6.4) 


where |i) is the initial state and |/) the hnal state after the absorbing a gluon. For our 
problem with a mixed-color heavy quarkonium system in the QGP, we assume that A ^ A 
in Eq. (|2.4|) and so the initial and hnal states are characterized by Oq = 1/\/2 = —as- Then 
we obtain 

47ragQ = (6.5) 


For these mixed-color states, if we represent the interaction between Q and Q by a gluon 
exchange, the associated coupling constant will be 5 f^[(—4/3) -|- (l/6)]/2. We identify this 
coupling constant with the coupling constant extract from lattice calculations, |(—4/3)aefr| 
of Eq. (|3.2|) . Hence, the coupling constant for the absorption of a gluon is then 


1 12 

ttoO = — X -— X 

sy 6 7 


4 

—-tteff 


( 6 . 6 ) 


Previously, the cross section for the dissociation of a color-singlet heavy-quarkonium in 
collision with a gluon was obtained in the short-distance approach by Peskin and Bhanot 
P5| , |2B[] and subsequently used by Kharzeev and Satz [^, for the simple case when the heavy- 
quarkonium can be treated as a hydrogenic system. The calculation of Peskin and Bhanot 
uses the operator product expansion, dispersion relations, and the sum of a large number 
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of diagrams in the large Nc limit. The leading order cross section involves the square of the 
color electric held which corresponds to the dissociation by an electric dipole radiation. 
As the physical process included in the evaluation of the dissociation cross section in Peskin 
et al. p5|,^ and in Eq. (|6.1|) of Blatt and Weisskopf are the same, they should give the 
same dissociation cross section. We shall show the equality of both results explicitly for the 
simple case considered by Peskin and Bhanot in the Appendix. 

The analytical result of Peskin and Bhanot is useful only for the idealized case of a 
hydrogenic 15 initial color-singlet state and a non-interacting IP hnal state in the large N^. 
limit. The wave functions of the heavy quarkonium systems of interest are not hydrogenic. 
The energy and the wave functions of the initial and hnal states also depend sensitively on 
temperature. The large limit is only a crude approximation for = 3. We need to use 
the more general results of Eq. ( p.l|) for our problem in which both the initial and the hnal 
state wave functions have been obtained numerically, without using the large Nf. limit. 



Fig. 6. The El cross section for the dissociation process T(15) -|- 
g ^ h + h after absorbing a gluon. 

With the L = 1 phase shift 5i(r) and amplitude ai(r) obtained in Section IV and the 
L = 0 bound state wave function obtained in Section III, the overlap integral (|6.2|) can be 
evaluated to give the dissociation cross section of T(15) by gluon absorption. We show in 
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Fig. 6 the dissociation cross section of T(IS') by gluon absorption as a function of the hnal 
continuum state energy e. The cross section is zero at e = 0, rises to a maximum and then 
decreases with increasing e. The dependence of the El T(IS') dissociation cross section on 
energy is similar to that of the El deuteron photodissociation cross section as shown in Fig. 
XII.4.1 of El. 


Because the mass of the b quark is large, the Ml cross section is small relative to the 
El cross section, and can be neglected. 

It should be noted that the dissociation cross section has been calculated for the exclusive 
channel g + T{lS) —^ b + b. As the energy increases, the cross section for this specihc channel 
decreases (Fig. 6) and other reaction channels such as g + T{lS) —> b + b + ng + n'{q + q) will 
make important contributions to the dissociation of the T(IS'). It will be of great interest 
to evaluate these cross sections in future work. 


VII. DISCUSSION 


We have studied the dissociation of a heavy quarkonium state in a quark-gluon plasma 
using the hnite-temperature potential VQQ{r) inferred from lattice gauge calculations of 
Karsch et al. Focusing our attention on states that can be directly or indirectly detected 
by dilepton decay, we can classify quarkonia into three groups: (1) quarkonia which dis¬ 
sociate spontaneously in hadronic matter: -0', Xci Xbi (2) quarkonia which 

are stable in hadronic matter but dissociate spontaneously in the quark-gluon plasma: Xb'-, 
and (3) quarkonia which are stable in hadronic matter and the quark-gluon plasma below a 
dissociation temperature: T(IS') with the dissociation temperature of 1.1 IT^. 

From this classihcation, one notes that the Xbi and Xb 2 mesons are good indicators of a 
quark-gluon plasma, as they are stable in hadronic matter but dissociate spontaneously in 
a quark-gluon plasma. The T(IS') (and p;,(lS')) have the distinction that they are the only 
stable quarkonia in the quark-gluon plasma, although the binding energies are rather small. 

For a quark-gluon plasma at a temperature below the T(IS') dissociation temperature 
Td, a direct observation of the position of the T(IS') using decays to dileptons is of great 
interest. The mass of the T(IS') state in the quark-gluon plasma is 2mb + e(l*S'), where 
e(lS') lies between -0.0175 GeV and 0 GeV (Fig. 3), and the current b quark mass mb is 
4.2 ± 0.2 GeV |^. The invariant mass of the dilepton peak from the decay of the T(IS') in 
the quark-gluon plasma will be substantially lower than the dilepton invariant mass of 9.46 
GeV seen in the decay of the T(IS') in free space. The fraction of dileptons from the decay 
of the T(IS') inside the plasma depends on the lifetime of the quark-gluon plasma and the 
dilepton decay partial width of T in the plasma. This fraction may be rather small, making 
the observation of the dilepton energy peak shift a difficult task. 

We have examined the dissociation cross section for the exclusive channel g + T(IS') —> 
b + h through a color El transition. We have also studied the phase-shifts for Q-Q collisions 
in the continuum. We find that there are T(IS') and rjbi^S) resonances at e = 0 for L = 0 
at the dissociation temperature Td = l.llTc, with a half width of F/2 ~ 0.2 GeV. The 
occurrence of these L = 0 resonances just above e = 0 is a characteristic of the quark- 
gluon plasma, as they are not expected from other sources. These resonances will decay 
in the quark-gluon plasma and appear as separated b and b quarks. After the quark-gluon 
plasma cools to become hadronic matter, the b and b quarks will hadronize into B and 
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B mesons. The invariant mass of this pair will show a resonance at low relative kinetic 
energies. However, direct detection of B and B mesons is a difficult task whose prospects 
may improve with future advances in detector technology. 

The author would like to thank Drs. M. Creutz, Keh-Fei Liu, L. McLerran, and S. Ohta 
for valuable discussions on lattice gauge theory. The author also wishes to thank Drs. T. 
Barnes, V. Cianciolo, H. Crater, and G. Young for helpful discussions. This research was 
supported by the Division of Nuclear Physics, Department of Energy and the Laboratory 
Directed Research and Development Program of Oak Ridge National Laboratory under 
Contract No. DE-AC05-00OR22725 managed by UT-Battelle, LLC. 

Appendix : Approximate gluon dissociation cross section for a color-singlet 
heavy quarkonium 

We would like to show how Eqs. lead to the result of Peskin and Bhanot 

P3| , |2B[] for the special case they considered. They assumed that the initial ground state 
of the color-singlet heavy quarkonium can be approximately described by a hydrogenic IS* 
wave function 


Mi5(r) = 27^/Ve (1) 

where 7 ^ = 2hqqB, B = /iQ<5asingiet/2, and 7 = l/iggOsingietl- Here, Osingiet is the color-singlet 
coupling constant, relating to the strong interaction coupling constant by 

- 1 

Q^singlet = <^s, (2) 

where = 3 is the number of colors. Peskin and Bhanot used the large A^ approximation 
and the above color-singlet effective coupling constant for A,, = 3 is approximated in this 
large limit by 

r -_3 

^singlet 

large Nc 2 2 


Peskin and Bhanot considered the hnal color-octet state to be approximated by the IP wave 
function as in the free-particle case. The corresponding wave function is 


Mip(r) 


sin kr 

—;-cos kr. 

kr 


The overlap integral I of Eq. ( 0 ) then gives 

pOO 

I = uip{r)uis{r)rdr 
Jo 

1675/2/^2 

= (p+y)3- 


( 4 ) 


( 5 ) 


Using /c^ = 2fiQQe and the relation between 7 and the binding energy P, Eq. ( |0|) from 
Blatt and Weisskopf gives the dissociation cross section 


'^d]^{iQQ)is + 9 ^ Q + Q) 


47^. .2 «gQ je/Bf/'^ 
3^ ’ + 


( 6 ) 
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The quantity OgQ in the above equation is the effective coupling constant for the interaction 
of a gluon with the Q (or Q) of the quarkonium. In the problem studied by Peskin et al. 
p5| , |26[] , the initial IS' state is in the color-singlet state 

= ( 7 ) 

V3 y 

and the final IP state is in a color-octet state with color component c 

|8c) = (8) 


The fermion-(gauge boson) vertex in a Feynman diagram is associated with —in QED 
and with g'y^\^/2 in QCD where c is the color index of the gluon. Thus, in going from 
the QED case of Blatt and Weisskopf to the QCD case of Peskin and Bhanot, e in QED is 
replaced by fi'A‘^/2 in QCD. The value of = dyragEo in QED of Blatt and Weisskopf is 
replaced in QCD by 


I(8c|c/y| 1)P = dvragQ 


where acts only on the Q (or Q). Eqs. gives 

®gQ ~ 


9 


dvr X 6 


Otg 

~6 


(9) 


( 10 ) 


The above equations lead to 


o'{{QQ)is + 9 ^ Q + Q) 




IbvrA 1 (e/P)^/^ 

3c/V (e/P + l)5’ 


( 11 ) 


which is the same as Eq. (4.4) of Bhanot and Peskin (Ref. [^), demonstrating the equiva¬ 
lence of the descriptions of Blatt and Weisskopf and Peskin and Bhanot. 

It should be pointed out that the large W limit is only a crude approximation for W = 3. 
If one does not invoke the large W limit, we then have for W = 3, Osingiet = —(4/3)as. Hence, 
we obtain a more accurate dissociation cross section from Eqs. (|6D-(pTD 


o'{{QQ)is + 9 ^ Q + Q) 



IGttA 1 (e/P)^/^ 

vj ^ (e/P + 1)^’ 


( 12 ) 


which is (9/8)^ times the approximate large-W limit result of Peskin and Bhanot 


One notes that from this cross section expression (0), the peak of the dissociation cross 
section is located at e = 3P/7 and the maximum dissociation cross section from Eq. ([T^) is 


cr„ 


i{iQQ)is + 9 ^ Q + Q) — 


18.90 


(13) 


If one uses the hydrogenic description for J/^jj with = 0.3, then P = 0.06 GeV and 
the maximum of the gluon dissociation cross section is located at e = 0.025 GeV with a 
maximum of 10.90 mb. For a color-singlet T with as = 0.2, then P = 0.089 GeV and the 
maximum cross section is located at 0.038 GeV with a maximum of 1.5 mb. 
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